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We analyse the accuracy of the approximate WKB quantization for the case of gen- 
eral one-dimensional quartic potential. In particular, we are interested in the validity of 
semiclassically predicted energy eigenvalues when approaching the limit E — > oo, and in 
the accuracy of low lying energy levels below the potential barrier in the case of generally 
asymmetric double-well quartic potential. In the latter case, using the standard WKB quan- 
tization an unnatural localization of eigenstates due to the negligence of tunneling is implied 
and thus the validity of semiclassics is uncertain. In all computations the higher order cor- 
rections to the leading semiclassical approximation are included using the complex contour 
integration technique. We show that these corrections can improve accuracy of semiclassical 
approximation greatly by many orders of magnitude. 

§1. Introduction 

The semiclassical approximation is a very useful approximative tool in quantum 
mechanics, not only because it often yields analytic expressions but also because it 
helps us to understand the correspondence between quantum features and classical 
mechanics, that is particularly important in the quantum chaos. 

Certain rather naive thoughts about semiclassical approximation have been es- 
tablished in the past. Namely, it has been believed that the semiclassical quanti- 
zation gives sufficiently accurate energy values of quantum eigenstates, when h is 
sufficiently small, in comparison to some action that stands with it, and in the case 
of high energies as well. Both limits imply that the de Broglie wavelength must 
be sufficiently small, but they are nevertheless not equivalent unless the underlying 
Hamiltonian has the scaling property. But the crucial criterion for assessing the 
accuracy of semiclassically predicted energy values is the comparison of the abso- 
lute value of the error of semiclassical energies with the mean energy level spacing 
(Prosen and Robnik 1 )). The fine structure, the statistical and other properties of 
energy spectra cannot be reliably resolved unless the error is small enough compared 
to the mean energy level spacing. 

The usual semiclassical techniques like EBK quantization (Maslov 2 ) ) and Gutz- 
willer trace formula are just the first (leading) terms of a certain semiclassical expan- 
sion of quantum quantities in the h~ power series. Typically these are only asymptotic 
series, but occasionally can nevertheless be convergent, as a surprise (see below). The 
absolute error of energies, predicted by those two methods is typically of order h 2 , 
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while the mean energy level spacing scales as ti* , where / is the number of degrees 
of freedom (dimension of the system's coordinate space). Thus, the two methods 
seem to be applicable only in systems whose dimensionality is less than three. At a 
fixed value of Ti and in the limit of high energies E — > oo, the situation is even worse. 
The error of the leading order semiclassical approximation, measured in the units of 
the mean energy level spacing, can diverge or is bounded from below as E goes to 
oo even in cases of one and two dimensional systems (Prosen and Robnik 1 ), Rob- 
nik and Salasnich 3 )' 4 )). Nevertheless, for some specific systems, all of them exactly 
solvable, the WKB series can be convergent and its sum equal to the exact result. 
Seeref. 3 )' 4 )- 7 )- 10 ). 

In this paper we analyse a one-dimensional general quartic potential, for which 
we are able to calculate the systematic higher order correction in the ft expansion. 
Firstly we examine the simple case of pure quartic potential, also called homoge- 
neous quartic oscillator, which is a scaling system, and thus we expect the fi ^ 
limit, to be equivalent to the E — ► oo limit, so that the error of semiclassical en- 
ergies in units of the mean energy level spacing should converge to zero as energy 
increases. Secondly, we tried to enrich our understanding of global behaviour of the 
semiclassical approximation by examining its validity in the case of potential with 
more than one minimum, like asymmetric double-well quartic potential. A stan- 
dard semiclassical quantization condition, obtained by requiring that the complex 
wavefunction as a function of the complexified coordinate be singlevalued, demands 
that the change of the semiclassical phase around a circuit enclosing the two turning 
points (at the boundaries of the classically allowed regions) is an integer multiple 
of 2ir. We assume that neglecting the tunneling, as an exponentially small effect, 
this quantization condition can be also used in the cases where classically allowed 
regions are separated in the coordinate space with the potential barrier. In these 
cases motion in each well is meant to be semiclassically quantized separately, and 
the error of energy values because of this unnatural localization is supposed to be 
exponentially small. Strictly speaking, we actually do not know rigorously how to 
interpret this quantization condition, which is another important reason to study the 
implementation of WKB methods and to test them in multi-well potentials like in 
our present work (also see Robnik and Romanovski 9 ) ) . However, if one is interested 
explicitly in the fine structure of energy spectra, especially and in particular the 
exponentially small splitting (due to tunneling) in symmetric double well potentials, 
then the approach expounded in Robnik et al 5 ^ is the appropriate one. 

We show that the homogeneous and asymmetric nonhomogeneous quartic po- 
tential is indeed an example of a potential where the semiclassical approximation 
works well in all energy ranges that we were interested in. By adding higher order 
corrections, the accuracy of semiclassically predicted energy values measured in units 
of the mean energy level spacing can be increased by many orders of magnitude. 
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§2. Homogeneous quartic potential and semiclassical methods 

The Hamiltonian of the simple homogeneous quartic potential, that we consider 
in this section is 

H = ^- + Ax 4 = -^^- 2+ Ax\ (24) 
2m 2mox z 

Introducing the new scaled and dimensionless variables 

1/3 

/ mji. \ -' ' ^ _ / rrr \ 

x 



(if) • e=e Wa) • R=1 - {M> 



the Hamiltonian can be rewritten in the simple form 

k = ~\^ +il - < 2 ' 3 > 

In all further calculations only the scaled form of the Hamiltonian and of the variables 
will be used, so that the tilde will be omitted in the further writing. Note that the 
scaled energy is the only one system parameter. 

We study this system for two main purposes. Namely, we are interested in the 
asymptotic (E — ► oo) behaviour of standard semiclassical quantization formulas of fi- 
nite order in the case of scaling system like this, and further we want to compare the 
spectral properties of homogeneous quartic potential with the spectral properties 
of general asymmetric double well quartic potential which we also call non-trivial 
potential. Again, in this comparison we are mainly interested in the asymptotic 
properties meaning high energy behaviour of the spectrum of a non-trivial potential 
whose geometry is asymptotically the same (as |x| — > oo) as for the homogeneous 
quartic potential. For these reasons we use the standard WKB quantization proce- 
dure rather than the implicit semiclassical formulas of Voros 6 ** . 

To calculate the semiclassical spectrum of a Hamiltonian we observe that the 
wavefunction can always be written in the form 



Y>(x) = exp (j:<t{x)) , (2-4) 



where a(x) is the complex phase that satisfies the Riccati differential equation (for 

a'\x) + V'(ar) = 2m(E - V (x)) . (2-5) 
The WKB approximation rests upon the expansion of the phase in the fi power series 

*(*) = E 7 (2-6) 



k=Q 



Substituting (2-6) into (2-5) and comparing like powers of fi gives the recursion 
relation for a' k (x) (n > 0) 



a' 2 = 2m(E-V(x)), £ aVn-fe + o" n -i = . (2-7) 

fc=0 
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The quantization condition is obtained by requiring that the complex wavefunc- 
tion (2-4) is singlevalued as a function of the complexified coordinate, which means 
that the change of phase, when going around the classically allowed region between 
the two turning points a (V(a) = E) and b (V(b) = E), should be equal to an integer 
multiple of 2irh 



k=0 



<j> da k = 2TThn. (2-8) 



The integer n is the quantum number. The zeroth and the first order terms give 
the EBK or torus quantization formula, where the first order term gives the Maslov 
corrections with Maslov index 2 for one dimensional system. We observe that all the 
other odd terms in (2-8) vanish when integrated along the closed contour around the 
pair of turning points, because they are exact differentials which is a highly plausible 
hypothesis yet to be proved rigorously (Bender et al 7 \ Robnik and Romanovski 9 ) ) 
but it has been demonstrated to hold for some low orders. We should stress that 
the closed contour is the contour in the complex x plane, since all integrals in (2-8), 
except the zeroth order term, diverge when integrated on an interval along real x- 
axes between the two turning points (Bender et al 7 \ Robnik et al^). However, there 
is a way of expressing these integrals as a partial derivative w.r.t. the energy of a 
fundamental real integral on the interval between the two turning points. Finally, 
taking this into account, the semiclassical quantization condition can be written as 

]T (jj j da 2k = 2nh (n + - J , (2-9) 

that is the sum over even-numbered terms only. The first three terms in (2-9) are 

j da = 2v / 2^ £ \Je~ V(x) dx, V(a) = V(b) = E, (2-10) 

d r b V"(x) 

^- / - { ' dx , (2-11) 

uV^dEJa ^E - V(x) 



dao 



do~4 = 



7 8 3 f b V" 2 (x) , 8 2 f b V""(x) , ' 
dx — - — / — , = dx 



576^2^3 [5 d E 3 J a y/E - V{x) ' dE 2 J a y/E-V{x) 

(2-12) 

A straightforward calculation of these terms for the case of the homogeneous 
quartic potential gives the semiclassical quantization formula with corrections up to 
the fourth order 

fc=2 



E 3 /* J2 a k E- 3k ' 2 = Jn + \ 
k=o v 1 



a- F ^ 2 a- ^ a- 1 -^^- (2 13) 

a °-^7^' 1 ~~^f } 2 "^7^' ( 3) 
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Fig. 1. The error of the 1100 lowest semiclassically predicted eigenenergies of the homogeneous 
quartic potential (HQP). Negative decadic logarithm of the error in units of mean energy level 
spacing (-Log dE) is plotted against energy E. The nearest level spacing is averaged over 7 
nearest states to obtain the mean energy level spacing. 



The exact spectrum is obtained by numerical diagonalization in QUADRUPLE 
PRECISION (32 valid digits) of the Hamiltonian matrix in the basis of eigenstates 
of harmonic oscillator, whose Hamiltonian is 

H = \{f + ^). 

The matrix elements in this basis are 

H n ,m = (^(2n + 1) + ^ (2n 2 + 2n + 1)) 5 n , m + 

+ ^/(n+l)(n + 2) (-j + -j^ ( 2n + 3 )) S n+2,m+ 

+^/(n + l)(n + 2)(n + 3)(n + 4)^6 n+4 , m . (2-14) 

The results of our calculations are plotted in the figure 1, where the negative 
decadic logarithm of the error of the semiclassical energies obtained from the zeroth 
order, the second order and the fourth order semiclassical quantization formula in 
units of the mean energy level spacing is plotted against energy. In table I of the Ap- 
pendix B we show the numerical values for some representative energy levels, showing 
the fast increase in accuracy with higher orders of the semiclassical approximation. 

^From the plots we can deduce that the error really goes towards zero as we go 
to the high quantum numbers and/or energies. The accuracy of the semiclassics can 
be vastly enhanced, when higher order corrections in fi expansion are added to the 
EBK or torus quantization. Contrary to the case of potential Uq/ cos 2 (ax) previously 
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studied (Prosen and Robnik 1 ), Robnik and Salasnich 3 )), the homogeneous quartic 
potential is the case, where the ft — > limit is equivalent to the E — > oo limit as we 
predicted in the introduction, due to the scaling properties of our system. 

§3. Asymmetric double-well quartic potential 



The Hamiltonian of the system that we consider in this section is generally 
written as 

° 2 + Ax 4 + Cx 3 -2Bx 2 . (3-1) 



H 



2m dx 2 

Again we introduce the scaled variables 




h 



CAB* '" eff ~"VmB-^ 
and rewrite the Hamiltonian containing these variables 

d 2 



h- 



A 



H 



h eff 



2 dx 2 



+ x 4 + Xx 3 - 2x 2 . 



(3-2) 



(3-3) 



The notation with tilde has already been omitted in the equation (3-3). The 
H e ff and A are two additional parameters of the asymmetric double- well quartic 
potential. The latter is called the parameter of asymmetry. To check the accu- 
racy of the semiclassical quantization for this type of the system, we perform the 
semiclassical quantization separately for states with energy below the potential bar- 
rier (E < 0), and for the states with energy above the potential barrier (E > 0). 
Neglecting the tunneling, those states 
below the potential barrier, are semi- 
classically quantized using the stan- 
dard semiclassical procedure, de- 
scribed in the previous section, for 
both wells separately. We assume 
that the error committed by this ap- 
proximation and localization is only 
exponentially small, since the tun- 
neling is only exponentially small. 
The calculation of the integrals (2T0), 
(2-11) and (2-12) that enter the semi- 
classical quantization formula (2-9) is 

much more complicated here than for the system discussed in the previous sec- 
tion. The expressions obtained immediately after the integration on coordinate x 
in (2-11) and (2-12) are very complex functions of the four roots of the momentum 
p, and energy. These functions should be further differentiated with respect to en- 
ergy E. Please note, that the roots of p = \/2(E — V(x)) are also dependent on 
energy and in case of the fourth order correction formula (2-12) a three-fold differ- 
entiation with respect to E should be done. The expression for the fourth order 
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Fig. 2. A sketch of the asymmetric double-well 
quartic potential (ADWP) 
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corrections is thus too complex to be handled with our hardware equipment, re- 
gardless of the fact that we do not need to know the explicit dependence on E for 
the roots. The derivative of a root a with respect to the energy could be written 
in a form da/dE = (dV(x)/dx)~ 1 \ x=a . In most of the manipulations we used the 
Mathematica software. The calculation of expressions (2T0) and (2T1) are briefly 
sketched in the appendix A. 

We have performed our calculations for three chosen values of h e ff, namely at 
h e ff = 0.01, h e ff = 0.1 and Ti e ff = 1 5 an d at a fixed value of the parameter of 
asymmetry A = 0.2. The shape of our asymmetric double well potential is shown in 
figure 2. The negative decadic logarithm of the error of the energies calculated with 
the zeroth order and the second order semiclassical quantization formula, in units of 
the mean energy level spacing, is plotted in the figures 3 and 4. 

As we see from figures 3, 4 and from the table II in the appendix B, the accuracy 
of semiclassical methods increases with energy like in the homogeneous case. The 
accuracy is about 4 orders in magnitude or more greater when the second order 
corrections are included in the semiclassical quantization formula. 

It is quite surprising how accurate the semiclassical quantization is even in the 
range of the energies below the potential barrier whilst close to the top of the barrier 
it is not so good, as expected. Of course we expect that the semiclassical methods are 
valid where the potential is similar to the harmonic potential and our expectation 
seems to be proved now. As we see the accuracy decreases and the improvement 
of the semiclasically predicted energies by considering the higher order semiclassical 
corrections e.g. a 2 gets smaller as we approach the top of the barrier. This is 
expected, since the exponentially small error caused by the negligence of tunneling 
increases with energy, as we have shown in our previous work (Robnik et al 5 ^). 
Our results of this section are depicted in figures 4-6. Generally, the accuracy of 
semiclassics increases as we decrease the value of Ti e ff, what is of course expected 
and shown in figure 6. 

Again, the exact spectrum has been calculated by numerical diagonalization in 
QUADRUPLE PRECISION of the Hamiltonian matrix in the basis of eigenstates 
of the harmonic oscillator. Because the symmetry w.r.t. parity operation is broken 
here, the matrix does not decompose into even and odd part any more. The matrix 
elements of the matrix in the same harmonic oscillator 




basis as before read as 




+3A 
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(n + l)(n + 2)(n + 3)<5 n+3 , m + 
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Fig. 3. Negative decadic logarithm of the error of the semiclassical energies in units of the mean 
energy level spacing (-Log dE) vs. energy (E) of 200 lowest states (left), of which 9 lie below the 
barrier and are separately plotted in the right graph. Asymmetric double- well quartic potential 
with h e ff =0.1 and A = 0.2. 




Fig. 4. Negative decadic logarithm of the error of the semiclassical energies in units of the mean 
energy level spacing (-Log dE) vs. energy (E) of 200 lowest states, of which 86 lie below the 
barrier. Asymmetric double- well quartic potential with h e ff = 0.01 and A = 0.2. 
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Fig. 5. The mean energy level spacing (DE), averaged over 7 nearest states, vs. energy (E) of 200 
lowest states. Asymmetric double- well quartic potential with h e ff —0.1 and A = 0.2 (left) and 
h ef f = 0.01 and A = 0.2 (right). 
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Fig. 6. Negative decadic logarithm of the absolute error of the semiclassical energies of the ground 
state (-Log d(aE)) vs. h e ff (h). Asymmetric double-well quartic potential with A = 0.2. 



§4. High energy correspondence between the exact spectra of the 
homogeneous and the general quartic potential 

During our work new questions about the spectral properties of the systems 
with non-trivial potential have arisen. By the word non-trivial potential we mean 
e.g. the asymmetric double-well quartic potential and similar potentials with compli- 
cated structure in certain regions of coordinate space and with a simple asymptotic 
behaviour as \x\ — > oo. The question is, whether the spectra of the homogeneous 
quartic potential and that of the asymmetric double-well potential with Ti = 1 con- 
verge pointwise or at least to a certain shift in energy and in the quantum number, 
as E — > oo. 

The scaled Hamiltonian of the general quartic potential can be written as 

# =y +X 4 + flV 1 (x) 



Vi(x) =Xx 3 -2x 2 . 



(4-1) 



Accuracy of the WKB approximation 



11 



In the following figures we present the comparison of the spectrum for the case 
ix = with the spectrum of the system where \i = 1, A = 0.2. We compare about 
500 lowest states. In the figures 7 the mean energy level spacing against energy is 
plotted for both systems. The discrepancy between the curves for the two systems 
is clearly seen. To any state of \i = case we find the state of \x = 1 case which is 
energetically most close to it and in figure 8 we plot the discrepancy of those energies 
measured in the units of the mean energy level spacing of the \i = case, against 
the energy of ji = states. We clearly see some cut-offs that occur at some energy 
values, where the difference in quantum numbers of the energetically nearest states 
of the two systems is enlarged by 1 . 

To understand these features we calculate the zeroth order semiclassical approx- 
imation for the mean energy level spacing DE (the Thomas- Fermi formula). 



DEadwp = -r~ = k 
an 



L 



Q2 dx 
<n p(x) 



, p=p(E-V{x)), p(oi)=p(a 2 ) = 0. 

(4-2) 

For the case of the asymmetric double well potential (ADWP) the integral in the 
expression above becomes 



L 



dx 



ai y/E-V(x) 



F(m) 



(4-3) 



k 2 = (ai - b) 2 + c 2 , q 2 = (a 2 - b) 2 + c 2 , m = 



kq - (ai - 6)(a 2 - b) 
2kq 



where b and ±c are the real and the imaginary part of the other two complex conju- 
gate roots of p(x) (see appendix A). To calculate the mean energy level spacing for 
the ADWP case we expand the 4 roots of E — V(x) into the energy E asymptotic 
series. The first few terms of those series are 



ai,2 = ^ - ^ T 



3(A/x) 2 ii 
32 2 



(A/x) 2 M 
32 4 



15(A/z) 4 + 160(A;/)V + 256;/ 2 3/4 / 5/4 x 
2048 ^ + °\ E ) ' 



(4-4) 



, , ■ Au 

b±ic= h Au 

4 



(Am) 2 /* 

32 4 



£-1/2 + i E l/4 



3(A//) 2 /x 



32 



+ 2 



f 15(A^) 4 + 160(A^) 2 /i + 256 j u 2 £ _ 3/4 
\ 2048 



+ o(V 5 / 4 ) 

Now, we can do the same with the quantities 2/y/kq and m 

32(A//)' 

Heq 128^/2 



£-1/4 



(4-5) 



2 = ^ 1/4 _ 3(A^) 4 + 32(A / ^ + 64 / ^_ 5/4 , 7/4 x 
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1 



m 



+ 



3(A//) S 
32 



+ 



/i 



+ o (V 3 / 2 ) 



Bearing in mind that 

<9F(m) 



3m 



2m(l 



m 



-(E(m) - (l-m)F(m)) 



we expand the elliptic integral F(m) into the asymptotic series around m 
Rewriting 

, , , r(i/4) 2 

F(l/2) - 



and 



E(l/2) 



8^ 



4^ 



1 + 



(4-7) 

(4-8) 
= 1/2. 

(4-9) 



8vr 2 



r(i/4) 4 



finally we get an asymptotic expression for the mean energy level spacing 



DE ADWP = E 1/4 j^ 



1 - E- 1 ' 2 



7T 



4T(l/4> 



: (16/z + 3(A/z) 2 ) 



(4-10) 



(4-11) 



It is easy to see that the factor outside the brackets is just the mean energy level 
spacing of the homogeneous quartic potential (HQP) 



DE HQP = E 1/A ^0^ 



(4-12) 



The zeroth order semiclassical predictions for the mean energy level spacing (4T1) 
and (4-12) are also plotted in figures 7. The result (4-11) can be obtained also 
perturbatively by a systematic perturbative treatment of the integral (4-2). 

Further we analyze the difference D of energetically nearest states of HQP and 
ADWP in units of mean energy level spacing of HQP. For this purpose we expand 
the zeroth order semiclassical quantization condition for ADWP 



ra-2 
/ 



pdx 



±V2 1 



{Xfi)na 2 k-a 1 q ( 3(A/x) 2 



k — q 



+ 



3Xfi 
~8~ 



16 



' (\f,) 2 \ (q 2 -Q 1 )(fc + g) 
»+— ^— U(n,m) 



f Ml 



n 



F(m)+ 

(4-13) 

(k - q) 2 
4kq 



(see the end of the appendix A) into the asymptotic energy series, valid at high 
energies. To do that we follow the similar procedure as in the case of derivation of 
the mean energy level spacing asymptotic expression (4-11). Here we also use the 
first order expansion of the three elliptic integrals for small arguments e and 5 



F(l/2 + e) = F(l/2) + (2E(l/2) - F(l/2))e , 
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E(l/2 + e) = E(l/2) + (E(l/2) - F(l/2))e , 
77(0 + <5, 1/2 + e) = F(l/2) + (2E(l/2) - F(l/2))e + 2(F(l/2) - E(l/2))<5 



in our case 



1 

m 

2 



3(A/x) 2 
32 



+ o (V 3 / 2 ) ; d = n = (V 3 / 2 ) 



and so we get an asymptotic expression of the condition (4-13) 



- | M + - 



3^ 



l + E- 1 / 2 



3vr 2 



4T(l/4)' 



-(16/x + 3(A//) 2 



+ 



(4-14) 



Note, that the leading term of the upper equation (4-14) is precisely the zeroth order 
semiclassical quantization condition for HQP (see (2-13)) 



7T 



£3/4 



3^/vr 



(4-15) 



Since the M-th energy eigenvalue E + AE of ADWP is energetically nearest to the 
A-th eigenenergy E of the HQP (by definition AE is smaller than one half of the 
mean energy level spacing), we expand the difference of equations (4-14) and (4-15) 
into the AE/E power series up to the first order and solve it for AE. Finally we can 
write an approximate expression for the energy difference of the energetically closest 
states of ADWP and HQP at high energies 



AE = DE AQWP 



(M -N) 



7T 



4r(i/4) : 



-(16^ + 3(A / u) 2 ) J E 1 / 4 



(4-16) 



of which the leading term in the units of mean energy level spacing of HQP is 



D = 



AE 



DE 



HQP 



(M - N) 



4T(l/4) 2 



(16/x + 3(A^) 2 )£ 1/4 . 



(4-17) 



The above expression is also plotted in the figure 8 together with the exact results. 

According to the expression (4-17) and also with respect to the exact results, 
plotted in figure 8 we can conclude that the two spectra do not converge pointwise 
and not even converge to a constant energy shift since the expression (4-17) as a 
function of energy never reaches the constant value of D at a fixed value of difference 
between the appropriate quantum numbers N and M of the homogeneous quartic 
potential and the asymmetric double well potential. 

The expression (4-11) that fits the exact data quite well indicates that the mean 
energy level spacing of the asymmetric double well potential converges to the mean 
energy level spacing of HQP as energy increases but it converges as slow as E~ l l 4 . 



§5. Conclusions 



According to our results we can conclude that the semiclassical approximation 
gives sufficiently good results in the case of the quartic potential. As we see the 
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Fig. 7. Mean energy level spacing DE vs. energy E of the states between 400th and 500th state of 
HQP and ADWP (p — 1, A = 0.2) (left). The right figure shows an enlarged window of the 
left figure at the energy range [5000, 5500]. Solid line is the zeroth order semiclassical prediction 
(Thomas- Fermi rule). 




1000 2000 3000 4000 5000 

E 



Fig. 8. Energy difference of energetically nearest states of fi = case (HQP) and fx = 1, A = 0.2 
case (ADWP) vs. energy E of about 500 lowest states of HQP. The difference or discrepancy 
D is measured in units of mean energy level spacing of the HQP. Solid line is the zeroth order 
semiclassical prediction. 

semiclassical energies in units of the mean energy level spacing converge to the exact 
ones for the case of homogeneous quartic potential as E — > oo. Due to this scaling 
properties of the system we can talk about only one limit in this particular system 
where the semiclassics becomes exact, namely the semiclassical limit which is equiv- 
alent to the h — > or E — ► oo limit. The accuracy of the semiclassical energies 
at the high energy values is not spoiled even in the case of asymmetric double-well 
quartic potential, which is not a scaling system. Even in the energy range below 
the potential barrier the accuracy of the semiclassics is very good, and it is bounded 
with the exponential error due to the negligence of tunneling. It is encouraging that 
the higher order corrections in the Ti expansion of the semiclassical quantization for- 
mula increase the accuracy by many orders of magnitude. Obviously, the corrections 
like these are needed in two or more dimensional systems to prevent the error to 
diverge or to be bounded from below as we approach h — > or E — ► oo limit. Some 
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recent related papers are Robnik and Salasnich 3 )' 4 \ Robnik and Romanovski 8 )' 9 \ 
Salasnich and Sattin 10 ) . Unfortunately we are not able yet to derive the systematic 
corrections of higher orders in fi expansion to the EBK quantization in more than 
one dimensional systems. Nevertheless this work helps us to understand the global 
behaviour and limitations of techniques based on semiclassical ideas. As the higher 
order corrections yield more accurate results in one dimensional systems we might 
hope that similar ideas will help to improve the accuracy of the semiclassics in many 
dimensional cases. 

In section 4 we discuss the asymptotic behaviour of exact spectra of nontrivial 
potentials. By the word nontrivial we mean that the potential in certain range 
of coordinate possesses a complicated structure, but asymptotically as \x\ — ► oo 
becomes similar to a simple (scaling) correspondent. We have a strong evidence, 
exact numerical and semiclassical, that the spectra of both systems do not converge 
pointwise to the same energy values and not even to a constant energy shift as the 
energy is increased. Even the convergence of the mean energy level spacing of the 
two systems is very slow when E — ► oo (as slow as E~ l / A ). 



Appendix A 

Calculation of the zeroth order and the second order terms in the 
semiclassical quantization formula 



In this section we show how to calculate integrals § dao and § da2 for the case of 
asymmetric double-well quartic potential. The integrals that enter the semiclassical 
quantization condition are calculated separately for the regions where the momentum 
p = \f%E — V(x)) has four and two real zeros. 

The contour integrals around the pairs of turning points entering the zeroth and 
the second order semiclassical formula (2-10) and (2-11) respectively can be written 
as 



Vxdx = - 



X" dx 



E 



dx 



12 



X 

1 
3 




dx 



x dx 



X = E- 

All integrals that we need are 

dx 



Xx i + 2x 2 



+ 



x 2 dx 



(A-l) 
(A-2) 



PI 
2 



X 



h 



x dx 



X 



x 2 dx 



and will be calculated separately for the case where X has four real roots and for 
the case where X has two real roots. 

A.l. Momentum p and X have four real roots (a < b < c < d) 

The integrals are written for the case of contour integration around a pair of 
turning points a and b. In case of integration around c and d only substitution 
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a < — ► c, b < — ► d has to be done. 

X = E - x 4 - Ax 3 + 2x 2 = (x - a)(b - x)(x - c)(x - d) 
Firstly, we perform a substitution 

c — dy 2 
x = ~\ T 

i - y 2 

in integrals Ii, I2 and ^3 and we get 

4 dy 4 



I 



\/{d-a){d-b) J yi h y 2 _ y 2^ y 2 _ y2 ) j (d - a){d - b) 



yi 



l b-c 
6-d' 



h 



y/{d-a)(d-b) 
-(d-c) [ 

Jv 



2/2 



1)2 



1 a — c 
a — gT 



2/1 ^(y 2 -2/1X2/2 -y 2 ) 



(i-y 2 )y/(y 2 -ylM-y 2 ) 

[dJi-(d-c)J 2 ], 



V(d-a)(d-6) 



/ 3 



;V2 



2/1 ^/(y 2 -yi)(y2-y 2 ) 
dy 



+ 



V(d-a)(d-6) 

/■J/2 

- 2d(d-c) / - _ 

tt-y 2 )y/(y 2 -ylM-y 2 ) 
Hd-cfT r^— " 

^ w (l-j/ 2 )Y(y2-j/?)(^-y 2 )_ 

= ^ ^ [d2jl " ^ " C)J2 + (d " • 

V(a - a)(a - 0) 

Factorizing the term 1/(1 — y 2 ) 2 , can be rewritten in the form 



J. 



(2) 



5* + i». 



dy 



dy 



rite 

J vi (1 -y) 2 ^{y 2 - yj) (y 2 -y 2 ) J m (1 + y) 2 : y \y 2 - y\){y\ - y 2 ) 
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To calculate S, we use the following relations 



w 2 = a y 4 + 4oiy 3 + 6a 2 ?/ 2 + 4a 3 y + a 4 



b { ? (y - c) 4 + 4b[ c > (y - c) 3 + 66 2 c; (y - cf + 46f (y - c) + i# 



u(c) 



(c), 



and 



(c) {y-cfdy (c) yw {y-c)dy 



& r 

J yi 



w 



2h 



(c) 



Jm 



dy 



lie) 
°4 



:'/2 



lyi (V ~ c)w " 4 J yi (y - c) 2 w ' 
w(yi) = w(y 2 ) = 0, c^(y 1 ,y 2 ) (A-3) 

(see Magnus a/ 11 )). 

Inserting = (y 2 — y 2 )(y 2 ~~ 2/ 2 ) an d evaluating the upper expression (A-3) for 
c = 1 and c = — 1 and subtracting what we get in both cases, we can express S as 



S = A 



b? - 26< 1 > 
2b? 



b\ 



(i) 



bY> 2b? 



^i 1 ) - f„? 



(i/i + 2/2) - vlvl - 1 ■ 



Integrals 



Jx = l F (m) , J 2 = — F(m) + — — , - 

yi y2 2/2(1 -yi) 



n(n,m), J 3 = y 2 E(m), 



yl-y\ 

1/2 



n 



1-y? 



can be found in (Gradshteyn ei a/ 12 )) or (Abramowitz and Stegun 13 )). Functions 
F(m), E(m) and n(k,m) are the complete elliptic integrals of the first, the second 
and the third kind as defined by Abramowitz and Stegun. 
Finally we can write 



da = 



8^2 



3 y/(d-b)(c-a) 



' 1 / 3A 2 \ , , „ 3cA / A 2 \ 

+ \ f 1 + lo") ( d -^)( c -«) E M+ 



F(m)+ 



+ t( 1 + y) ( c - & ) i7 (^ m )} ; 



(A-4) 
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(d — c)(b — a) b — a 
m = - — - — — , n = 



(c-a){d-b)' 



c — a 



and 



/ ' h2 = -T2WE ^i-l)(c- a) { i cd +ab + D F(m) " (d " b)ic ~ a) E(m) } ' 

(A-5) 

A. 2. Momentum p and X have two real roots (a\ < 02, b± ci) 



X = E - x 4 - Ax 3 + 2x 2 = (x - ai)(a 2 - x)({x - bf + (?) 



In this case the substitution 



x 



a>2 + a\ I tan 2 | 
l + ftan 2 | 5 



k 2 = (ai - b) 2 + c 2 



(03 - 6) 2 + c 2 



is performed in integrals I±, I 2 and ^3 and following the similar procedure as in the 
case of 4 real roots of X, the integrals can be rewritten in the following form 



F(m) ; 



m 



- (ai - 6)(a 2 -b) - c 2 
2kq 

a 2 fc — a\q 



n 



(k-q) 2 
ikq ' 



k — q 



„, x (a 2 -ai)A; + g 
F(m) 7, t n(n,m) 



akk — a 2 q„, , , „, , , , , fa 2 — ai ) fc + a „. 

-2- ^F(m + fcgE(m) - (ai + a 2 + 2b v -JT(n, m 

k — q 4: k — q 

By using these expressions, / aVo and j> d<7 2 can be finally written as 

8V2 1 



da 



g A a 2 /c — a\q f ^ 3A 2 ^ a\k — a\q 



4 k-q 



16 I k — q 



F(m)+ 



+ + kqE(m) + 

3A / A 2 \ (a 2 -a 1 )(k + q) „, ,1 



(A-6) 



A; 2 = (01 -6) 2 + c 2 , g 2 = (a 2 -6) 2 + c 2 , 
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m = 



kq — (ai — b)(a,2 — b) — c 2 



2kq 



n = — - 



( k-qf 
4kq 



and 



d(J2 = 



1 d 1 

~JidEjk~q 



a%k — a?o aok — a\q 2\ „. , , „. . 1 

(A-7) 



In the section 4 we use the generalized form of the asymmetric quartic potential 

2 

V(x) = x 4 + ^(Ax 3 - 2x 2 ) X = E - x 4 - Xfix 3 + 2fxx 2 = y . 
In this case the integral (A-l) is slightly modified to 



Xdx 



E 



dx 



Xfj, f xdx 



+ u+ 



3(V)_ 2 
16 



x dx 



(A-8) 



yielding the /x dependent expression (4-13) for do, which is identical to (A-6) for 
fi = 1. Please note that the parameters fc, g, m, n as well as the integrals Ii, I2 
and ^3 are expressed in terms of the four roots of X and momentum and so is their 
general dependence on \x already considered. 
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Appendix B 

Some tables of results 



Table I. Some eigenenergies of the homogeneous quartic potential. By n we denote the consecutive 
quantum number of the state, while Esxact is the exact energy value of the eigenstate, and 

E slmi' E Se m i> and E< slmi 

are the semiclassically predicted energy eigenvalues of the nth state, 
calculated with the zeroth order, second order and the fourth order semiclassical quantization 
formula. 



n 




771 

t^Exact 




E (o) 

Semi 




F (2) 
Semi 




E (4) 
Semi 








.6679862591 





.5462673250 





.6178440470 





.5994975436 






5577710827 




7809440985 




9454695175 




8186297455 


1 


2 


.3936440164 


2 


.3635614446 


2 


.4003574561 


2 


.3991273384 






8230311603 




0057710124 




2199197003 




3857621598 


2 


4 


.6967953868 


4 


.6705199316 


4 


.6968569596 


4 


.6965336443 






6364619622 




9480733314 




9157646544 




2176087907 


3 


7 


.3357299952 


7 


.3148026019 


7 


.3358836919 


7 


.3357508981 






2709279884 




0083816986 




7629989381 




3385710476 


4 


10 


.2443084554 


10 


.2265364336 


10 


.2443782077 


10 


.2443100572 






3877107597 




1101923867 




8948408161 




3952052722 


100 


643 


.1833913927 


643 


.1811392054 


643 


.1833914100 


643 


.1833913927 






8056430848 




7609028834 




9319827383 




8185826816 


101 


651 


.7305795719 


651 


.7283422013 


651 


.7305795888 


651 


.7305795719 






7519147803 




2623003831 




3668259301 




7642701540 


103 


668 


.9091219363 


668 


.9069134812 


668 


.9091219523 


668 


.9091219363 






1714901544 




9348071280 




2366319599 




1827702401 


105 


686 


.1986791466 


686 


.1964986907 


686 


.1986791618 


686 


.1986791466 






2251646322 




3792741559 




3254876433 




2354809723 


495 


5397 


.2204747691 


5397. 


.2196972873 


5397 


.2204747693 


5397. 


,2204747691 






2444635012 




4414211246 




7029195001 




2444710621 


496 


5411 


.7486521189 


5411. 


.7478756815 


5411 


.7486521192 


5411. 


.7486521189 






9592537080 




1512620174 




4045276381 




9592611981 


497 


5426 


.2865864811 


5426. 


.2858110844 


5426 


.2865864813 


5426. 


.2865864811 






2223840485 




4128681921 




6545729016 




2223914687 


500 


5469 


.9588010681 


5469 


.9580287730 


5469 


.9588010683 


5469. 


.9588010681 






5121968046 




5845185750 




9057033639 




5122040195 


1000 


13774 


.2520020553 


13774. 


.2515153776 


13774 


.2520020553 


13774. 


.2520020553 






0377443322 




5917115936 




4151998141 




0377446170 


1005 


13866 


.1108150810 


13866 


.1103300181 


13866 


.1108150811 


13866 


,1108150810 






7240132753 




5261252080 




0964842617 




7240135535 


1008 


13921 


.2992811341 


13921 


.2987970336 


13921 


.2992811341 


13921 


,2992811341 






4444693588 




5119763791 




8139930100 




4444696331 


1009 


13939 


.7076086294 


13939 


.7071248487 


13939 


.7076086295 


13939 


,7076086294 






9093649576 




4694282718 




2779132907 




9093652307 


1010 


13958 


.1220154946 


13958. 


.1215320331 


13958 


.1220154946 


13958 


,1220154946 






3736333880 




1516055405 




7412099391 




3736336599 
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Table II. Some eigenenergies of the asymmetric double-well quartic potential with h e ff = 1, A = 
0.2, (/i = 1). By n we denote the consecutive quantum number of the state, while EExact are 
the exact energy value of the nth eigenstate and E^ mi are the semiclassically predicted energy 
eigenvalues of the nth state, calculated with the second order semiclassical quantization formula. 



n 


EExact 


ir( 2 ) 

Semi 





-0. 17296458264338966970 


-0.1733140474651558 


1 


0.65867088851932470219 


0.7915823554184969 


2 


2.51962151763406920996 


2.501600214006679 


3 


4.63865676589239143594 


4.640576944608262 


100 


619.636315269976695643 


619.6363152899056 


101 


628.028868586785663262 


628.0288686061785 


103 


644.899655331227155754 


644.8996553496056 


105 


661.883432587416475526 


661.8834326048518 


495 


5329.38195574207411204 


5329.381955742333 


496 


5343.81915399412485822 


5343.819153994382 


497 


5358.26617031658854323 


5358.266170316844 


500 


5401.66599550160242755 


5401.665995501854 
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